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We theoretically consider the substrate-induced Majorana localization length renormalization in 
nanowires in contact with a bulk superconductor in the strong tunnel-coupled regime, showing ex¬ 
plicitly that this renormalization depends strongly on the transverse size of the one-dimensional 
nanowires. For metallic (e.g. Fe on Pb) or semiconducting (e.g. InSb on Nb) nanowires, the renor¬ 
malization effect is found to be very strong and weak respectively because the transverse confinement 
size in the two situations happens to be 0.5 nm (metallic nanowire) and 20 nm (semiconducting 
nanowire). Thus, the Majorana localization length could be very short (long) for metallic (semi¬ 
conducting) nanowires even for the same values of all other parameters (except for the transverse 
wire size). We also show that any tunneling conductance measurements in such nanowires, carried 
out at temperatures and/or energy resolutions comparable to the induced superconducting energy 
gap, cannot distinguish between the existence of the Majorana modes or ordinary subgap fermionic 
states since both produce very similar broad and weak peaks in the subgap tunneling conductance 
independent of the localization length involved. Only low temperature (and high resolution) tunnel¬ 
ing measurements manifesting sharp zero bias peaks can be considered to be signatures of Majorana 
modes in topological nanowires. 


I. INTRODUCTION 

Majorana fermions (MF), which were proposed the¬ 
oretically 80 years ago as real solutions of the Dirac 
equation in the context of understanding neutrinoJ^^, 
have recently found their incarnations in solid state 
system^2HIl as zero-energy localized excitations in topo¬ 
logical superconductors (TS). In addition to the defin¬ 
ing property of being their own anti-particles, MFs 
in solid state systems are known to have exotic non- 
Abelian braiding statistic^, making them of high the¬ 
oretical interest with a pote ntial application in topo¬ 
logical quantum computatiorP^I While the initial 
proposals invo lvin g exotic materials (e.g. p-wave 
superconductor^^, 5 / 2 -fractional quantum Hall stated, 
and topological insulatorPl) have so far escaped ex¬ 
perimental realization, more recent prop osals utiliz¬ 
ing spin-orbit-coupled semiconductord^l^ltlSl or magnetic 
adatomJSltSl proximity-coupled to s-wave superconduc¬ 
tors have been the major focus of experimental efforts, 
with several published works claiming the observation 
of signatures of MFs in such systemJi^H^. The ini¬ 
tial claim of the observation of “Signatures of Majo¬ 
rana Fermions in Hybrid Superconductor-Semiconductor 
Nanowire Devices” by Mourik et aZ.^ in the InSb/Nb 
system, following precise theoretical prediction^SHini^ -^as 
later exper iment ally replicat ed by several groups in both 
InSb/NtPsHUmi and InAs/hybrid structures, giv¬ 
ing considerable confidence in the universal nature of the 
underlying physical phenomena. Despite this apparent 
experimental success, the identification of the zero-bias 
peak as a MF is still debatecP^ more definitive proof, 
such as signatures of the predicted non-Abelian braiding 
properties, is yet to be detected as una mbig uous evidence 
for the existence of the Majorana mode^. 

The motivation underlying the semiconductor- 


superconductor heterostructure Majorana platfornP is 
the artificial creation of a spinless low-dimensional (ei¬ 
ther 2EISI11 or irjSHISl) p-wave superconductor supporting 
MF^. The effectively spinless p-wave superconductivity 
residing in the semiconductor serves as the TS here 
arising from a combination of spin-orbit coupling, 
spin-splitting, and ordinary s-wave superconductivity. 
The combination of spin-splitting and spin-orbit cou¬ 
pling in the semiconductor allows, under appropriate 
conditions (of large enough spin-splitting and spin-orbit 
coupling), ordinary s-wave singlet Cooper pairs to 
tunnel from the superconductor to the semiconductor 
enabling topological p-wave superconductivity with 
triplet superconducting correlation^^ to develop in 
the semiconductor through proximity coupling. The 
experi menta lly relev ant t opological system has been an 
InSiphiZliai or InA^lSmi nanowire on a Nb or A1 su¬ 
perconducting substrate with the Zeeman spin splitting 
achieved through the application of an external magnetic 
field. For the magnetic field larger than a critical value, 
which is given simply by the proximity-induced super¬ 
conducting gap in the nanowire (assuming the chemical 
potential in the nanowire can be taken to be zero), the 
nanowire becomes an effe ctiv e TS with zero-energy MFs 
localized at the wire end^^HHini These zero-energy MFs 
should lead to zero-bias conductance peaks (ZBCP) 
in the tunneling conduc tance measuremeniP, and the 
experimental observatioiP^E^lt^ of such field-induced 
ZBCP has been taken as evidence for the existence of 
MFs in these topological nanowires. 

Typically, these MFs are localized near the ends of the 
wire with a finite Majorana localization length ^ which 
equals the superconducting coherence length in the topo¬ 
logical nanowire. When the wire length L ^ the two 
MFs at the two ends of the wire are considered to be in 
the topologically exponential protection regime with the 
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MF wavefunction falling off as along the wire (mod¬ 

ulo some oscillations not of particular interest here^^HZSl) . 
The MF is a well-defined zero-energy non-Abelian mode 
only in this exponentially protected (L ^ regime. By 
contrast, for short wires or long coherence length (i.e. 
L < ^), the two end MFs hybridize, producing split peaks 
shifted away from zero energy (and these split peaks rep¬ 
resent fermionic subgap states rather than MFs), and 
the topological protection no longer applies. It is only 
when the Majorana splitting is exponentially small (i.e. 
L ^ the nanowire system can be considered to be 
topologicaP. Thus, the quantitative magnitude of ^ (or 
more precisely the dimensionless length ratio L/^) is a 
key ingredient in the physics of MFs. We note that the lo¬ 
calized zero-energy MF bound states are also often called 
the Majorana modes, and we use the terminology Majo¬ 
rana fermions and Majorana modes interchangeably in 
this paper to mean the same zero-energy MF subgap lo¬ 
calized TS excitations in a spinless p-wave superconduc¬ 
tor. The current experimental MF search is mostly fo¬ 
cused on looking for subgap zero-bias tunneling conduc¬ 
tance peaks associated with these localized zero-energy 
Majorana modes - ideally, the ZBCP should have the 
quantized value of 2e^//i, but experimentally the ZBCPs 
observed so far have actual conductance v alues fac tors of 
10 (10^) lower in semiconductor nanowire^^^E^E^ (ferro¬ 
magnetic nanowire^^. Finite temperature, finite wire 
length, finite tunnel barrier, finite experimental resolu¬ 
tion, unwanted fermionic subgap states, and possible in¬ 
elastic processes conspire together to suppress the ex¬ 
perimental ZBCP strength, and this non-observation of 
perfect ZBCP quantization, which is much more severe 
in the ferromagnetic nanowires than in semiconductor 
systems, remains an open question in the subject. 

The MF localization length (or equivalently the TS 
coherence length) ^ is often assumed to be given by 
the standard superconducting coherence length formula, 
^ ~ uf/A, where vp and A are respectively the Fermi ve¬ 
locity and the induced TS gap in the nanowire. This su¬ 
perconducting coherence length formula is certainly ap¬ 
propriate for MF localization if the nanowire can be con¬ 
sidered an isolated spinless p-wave superconductor with 
the Majorana modes localized at the wire ends (which 
serve as the defects localizing the MFs). But in the ex¬ 
perimentally relevant situation the nanowire is not iso¬ 
lated, it is in fact in contact with (or lying on top of) 
a superconducting substrate which provides the neces¬ 
sary proximity effect to produce the topological system in 
conjunction with spin-orbit coupling and spin splitting. 
The question therefore arises whether the MF localiza¬ 
tion length formula is modified from the simple coherence 
length formula, or equivalently, whether the Fermi ve¬ 
locity and/or the appropriate nanowire superconducting 
gap are renormalized by the substrate. This issue was in 
fact discussed by Sau et alW^ and Stanescu et al^ some 
years ago in the context of 2D sandwich structures in¬ 
volving semiconductor/superconductor and topological- 
insulator/superconductor heterostructures, and very re¬ 


cently by Peng et aZ.l^in the context of ID ferromagnetic 
nanowire on superconductor hybrid structures used in the 
recent Princeton STM experimentP^. (The actual system 
theoretically considered by Peng et al^ is in fact a he¬ 
lical magnetic adatom chain, not a ferromagnetic chain, 
on a superconducting substrate.) In the first part of the 
current work, we theoretically study the MF localization 
question in depth for ID topological nanowire hybrid sys¬ 
tems, discussing the substrate-induced MF renormaliza¬ 
tion for both semiconductor and ferromagnetic nanowires 
on an equal footing, comparing and contrasting the two 
situations. In particular, we address the important issue 
of how it might be possible that Majorana localization 
length could be strongly (weakly) renormalized in ferro¬ 
magnetic (semiconductor) nanowire systems studied in 
different laboratories. 

Recently, it has been proposed that the three ingre¬ 
dients for the semiconductor nanowire proposal, i.e. su¬ 
perconductivity, magnetization, and spin-orbit coupling, 
can be realized in ferromagnetic nanowires deposited on 
a spin-orbit coupled superconductor. Experimental ev¬ 
idence in the form of a weak and broad zero-bias peak 
seems to provide some sup port to this hypothesis. Sev¬ 
eral theoretical calculationS^^^lESI have shown that as a 
matter of principle Majorana modes can emerge in fer¬ 
romagnetic wires in superconductors, as had be en sug¬ 
gested in more mesoscopic geom etrie d^^ l ^^ l ^^ -*^. Moti¬ 
vated by earlier ST M workS*^, some of the theoret¬ 
ical workS^^UMMI have modeled the system to be a 
chain of magnetic atoms [i.e. Yu-Shiba-Rusinov (YSR) 
impurities] on the superconductor surface, with no di¬ 
rect hopping between the impurity orbitals. This class 
of propos als su pports MFs only in a limited parame¬ 
ter regimd^^RSI. Although the YSR limit and the fer¬ 
romagnetic nanowire limit are the two extreme crossover 
regimes (tuned by very weak and very strong inter-site 
hopping in the nanowire, respectively) of the same under¬ 
lying Hamiltonian (i.e. there is no quantum phase tran¬ 
sition separating them, it is simply a hopping-induced 
crossover from the YSR regime to the ferromagnetic wire 
regime as hopping increases), the TS properties in the 
two limits are very different. In the YSR limit, con¬ 
siderable fine-tuning of the chemical poten tial is neces¬ 
sary in order to achieve TS and MP32Eni^ whereas the 
TS with localized MF arise generically without any fine- 
tuning in the ferromagnetic nanowire limit of strong hop¬ 
ping. Thus, any possible generic existence of MF in the 
magnetic adatom chain on superconducting substr ates is 
more natural in the ferromagnetic wire limilPSEI] rather 
than in the YSR limit. Therefore, considering the sys¬ 
tem as a ferromagnetic wire in proximity to a supercon¬ 
ductor is the natural way to understand the zero-bias 
conductance. This puts the ferromagnetic chain and the 
semiconductor nanowire topological systems on an equal 
footing with the only difference being that in the semi¬ 
conductor wire (the ferromagnetic chain) case the spin¬ 
splitting arises from an externally applied magnetic field 
(an intrinsic ferromagnetic exchange splitting). However, 
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in the semiconductor nanowire proposal^^^ the decay 
length of the Majorana is typically found to be compara¬ 
ble to the bulk coherence length in the superconductor. 
This also seems to be a common feature in the simula¬ 
tions of the ferromagnetic nanowire systems so far since 
the substrate-induced renormalization of the nanowire 
parameters is not included in the theory, thus co nsid- 
ering the nanowire to be effectively isolatecP2l22lII]_ Qj^ 
the other hand, it was noted that YSR bound states in 
STM experiments appeared to show a m uch shorter de¬ 
cay length than the coherence lengtfpSl^II. Based on this, 
it was conjectured^ that the Majorana modes might ap¬ 
pear to be confined to length-scales shorter than the co¬ 
herence length because of the delocalization of the wave- 
function into the bulk superconductor. Very recent work 
using a spin-helical adatom chain modeP^has shown how 
this substrate-induced renormalization mechanism may 
suppress the coherence length, possibly leading to a short 
MF localization length in a magnetic chain on a super¬ 
conductor which is very strongly tunnel-coupled to the 
magnetic chain, even if the topological superconducting 
gap is very small. Anomalously short Majorana localiza¬ 
tion lengths were numerically demonstrated for similar 
models in Refs.l^ancP^. Whether such a scenario applies 
to the actual experimental situation of Ref.^^ is currently 
unknown. 

The actual MF localization length question is of great 
importance to the experimental observations ir(^ since 
the estimated TS gap in the Fe adatom chains on Pb sub¬ 
strates studied therein is very small 0.1 meV) leading 
to a rather long coherence length (or equivalently MF 
localization length) of > 100 nm (assuming no substrate- 
induced renormalization) which would be much larger 
than the typical length of the adatom chains (5 — 50 nm) 
used in Ref.l^. In such a situation, the TS system is not 
in the exponentially protected regime at all, and the two 
end MFs should hybridize strongly leading to ordinary 
uninteresting fermionic states at high energies. Thus, to 
the extent the observations iiP manifest MF features, 
one must understand how the very long TS coherence 
length (i.e. the MF localization length) associated with 
the small induced superconducting gap can be consistent 
with the existence of isolated (i.e. non-hybridized) MFs 
in a system where the wire length is shorter than the 
localization length. (This conceptually problematic sit¬ 
uation does not seem to arise in the semiconductor TS 
systemJi^l^ since the wire length (> 1 fim) is typically 
much larger than the MF localization length (^ 10 nm) 
in the semiconductor nanowire systems— in fact, system¬ 
atic experimental efforts appeared to have observed the 
predicted Majorana hybridization effect in the semicon¬ 
ductor nanowires in the regime of long MF lo calizat ion 
length induced by a large external magnetic fielcP^Eni ) jf 
the observed subgap conductance features in Ref.l^ are 
indeed implying the existence of MFs in the underlying 
ferromagnetic nanowire, as has been concludecP^, then 
it is imperative that a clear theoretical understanding is 
developed for why the TS coherence length being much 


larger than the topological wire length is not a prob¬ 
lem. One possible way out of this quandary, suggested 
by Peng et alW^ using a helical magnetic chain model, 
is a strong suppression of the MF localization length by 
the substrate so that the L 3> C condition is still satisfied 
because the TS coherence length is renormalized by the 
substrate. This solution comes with the caveat that that 
the strongly localized MF is accompanied by a power-law 
taiP^, which extends to the coherence length of the bulk 
superconductor. So while the localization length appears 
to be reduced on length-scales where the splitting is mea¬ 
surable in tunneling, it is not clear that this suppression 
significantly aids the splitting problem for quantum in¬ 
formation purposes because of the power-law tail. We 
emphasize that it is not clear at all at this stage that 
the condition (in particular, very strong tunnel coupling 
between the substrate and the chain) for such strong 
substrate-induced renormalization of the MF localization 
length is actually operational in the experiment of Ref.l^, 
but the possibility of the substrate-induced suppression 
of the MF localization length must be taken seriously 
since it provides a way forward for future experiments to 
test this hypothesis as a possible resolution of the long 
TS coherence length conundrum, e.g., by studying the 
MF hybridization or splitting systematically as a func¬ 
tion of the dimensionless ratio by changing L in a 
controlled manner. 

The possibility that the MF localization length is 
strongly suppressed by the substrate immediately brings 
up the question of whether such a substrate-induced 
renormalization phenomenon is also operational in the 
semiconductor TS systems and, if so, the possible im¬ 
plications for the semiconductor MF experiment^i^fl^fll 
which have so far been simply interpreted on the basis 
of the standard ^ ^ vp/A formula with no substrate- 
induced coherence length suppression. The corre¬ 
sponding renormalization question for semiconductor- 
superconductor^ hybrid structures was studied in 
depth in Refs.l^aiSS], }jgj.g .^g generalize the theory 
to ID semiconductors and ferromagnetic metals in prox¬ 
imity to bulk superconductors. One possible reconcilia¬ 
tion for why the ferromagnetic (semiconductor) nanowire 
MF localization length is strongly (not) renormalized 
by the substrate is simply by assuming that the ferro¬ 
magnetic adatoms (semiconductor wire) are (are not) 
strongly tunnel-coupled to the substrate superconductor, 
but we want to avoid such ad hoc assumptions. One ques¬ 
tion we address in the first part of the paper is how differ¬ 
ent the substrate-induced Majorana localization length 
renormalization can be in metallic and semiconductor TS 
nanowires assuming essentially identical conditions to be 
prevailing for the substrate properties (including equiva¬ 
lently strong tunnel coupling of the wire to the substrate) 
in both cases. 

In this paper, we discuss in general the localization 
length of Majorana modes in proximity-induced super¬ 
conductors. To set a context, we start by discussing 
the localization length of Majorana modes in the Kitaev 
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chairP in various parameter regimes finding that depend¬ 
ing on parameter values the localization length can vary 
from being of the order of a lattice spacing to more than 
many hundreds of lattice spacings as is typical for the 
coherence length in ordinary superconductors. Realiza¬ 
tions of such Kitaev chains where the Majorana local¬ 
ization length is of the order of several lattice sites have 
been proposed in quantum dot arrayP^^, and therefore, 
in principle, the variation in the MF localization prop¬ 
erties can be tested in the linear quantum dot arrays 
by suitably tuning the dot parameters. Following this, 
we explicitly consider the proximity effect of the bulk 
superconductor with a goal to understanding the length- 
scale problem for ferromagnetic nanowires in proximity 
to spin-orbit coupled superconductors as used in the ex¬ 
periment of Ref.^. While the superconducting proximity 
effect is simply modeled by a pairing potential in most pa¬ 
pers on the subject, a more microscopic consideratiorP^ 
suggests that the proximity effect should be represented 
by a non-local frequency-dependent self-energy 

S (w; r, r') = T (r) Gsc (w; r, r') T {r ')'^, (1) 


where the matrix T represents the hopping between the 
nanowire and the superconductor, and is 

the Green function of the superconductor. By approxi¬ 
mating the Green function by that of a bulk s-wave su¬ 
perconductor, we show below that the local part of the 
self-energy has the form of 

S(w;r,r)--^To + rT,^, (2) 

^sc 


where F is the parameter determining the strength of the 
superconducting proximity coupling in the nanowire. As 
we argue in Sec. |III E based on a microscopic deriva¬ 
tion, the proximity parameter F ^ {R/a)~^Ep, where 
a ~ 0.5 nm and Ep 1 eV are respectively a length 
scale on the order of the lattice constant and the Fermi 
energy in the superconductor, and R is the radius of the 
nanowire. In mesoscopic semiconductor nanowire geome¬ 
tries R ~ 20nm leading to F ~ O.OlmeV <C Asc, and 
this fits into the simple picture for the proximity effect 
where retardation effects associated with the frequency 
dependence in Eq. ([^ may be ignored. Atomistic fer¬ 
romagnetic wires^H are qualitatively different since R ^ 
0.5nm for these wires and the estimated F ^ leV. This 
clearly puts the analysis in the regime F ^ Asp, which is 
the strongly renormalized limilP^. Establishing this key 
difference between the MEs in semiconductor and metal¬ 
lic nanowires (i.e. the MF localization length is strongly 
renormalized in one, but not in the other, due to sub¬ 
strate renormalization arising from retardation effects in 
the proximity self-energy function) is a main goal of this 
paper. While it might appear that a proximity effect of 
F ^ Asc in the ferromagnetic nanowire system would 
produce an s-wave pairing in the nanowire that is much 
greater than l^sc in contradiction with experiment, the 
frequency dependence of the self-energy, where F enters 


as a parameter, ensures that the s-wave pairing in this 
case is Ag ^ ^sc, which is much less than F. 

Therefore, it is clear that the full frequency-dependent 
self-energy is critical to get the physics correct!^. 

We present calculations of the local density of states for 
the ferromagnetic nanowire in proximity to a spin-orbit 
coupled superconductor. We find that the frequency de¬ 
pendence introduces renormalization of all microscopic 
parameters in a way which drastically reduces the coher¬ 
ence length in the ferromagnetic nanowire, in contrast to 
the semiconductor nanowire MFs where the standard def¬ 
inition for the coherence length applies with little renor¬ 
malization by the substrate. After establishing that the 
MF localization in the ferromagnetic nanowire system 
could indeed be very short in spite of the induced topo¬ 
logical superconducting gap being very small, we consider 
the actual experimental situatiorPIl where the MF obser¬ 
vation in a Fe chain on a superconducting Pb substrate 
has been claimed in an STM study. At first, the very 
small induced gap (~ 100 ^eV) in the experiment seems 
to indicate a very long MF localization length much 
larger than the length of the Fe chains, casting serious 
doubt on the experimental interpretation of lattice-scale 
MF localization in the system. Peng et alW^, however, 
provided a way out of this puzzle by showing that, as¬ 
suming the Pb-Fe tunnel coupling is strong, it is possible 
for the MF localization length to be short in spite of the 
induced gap being small within their helical chain model. 
In the current work, we go beyond the analysis in Ref.^^ 
by showing that this renormalization is a generic feature 
of chain TS platforms. Moreover, we directly address 
the question of why such effects are weak in the semicon¬ 
ductor nanowire, revealing the key role of mesoscopic vs 
microscopic geometry in determining the substrate cou¬ 
pling strength. This provides a clearer understanding of 
the relationship between these the superficially different 
chain and nanowire TS platforms. 

In the second part of the paper, we address the issue of 
energy localization of the Majorana zero mode in the re¬ 
cent ferromagnetic wire STM experiment!^. Apart from 
spatial localization, another peculiar feature of the Majo¬ 
rana modes in the ferromagnetic nanowire experiment re¬ 
lates to the broadening of the spectrum in the energy do¬ 
main. This feature of the experiment!^ results from the 
high temperatures (T = 1.4 K) at which the experiment 
is carried out comparable to the induced gap 0.1 meV) 
itself. Such high temperatures limit the energy resolu¬ 
tion of the STM, thus compromising the claim of the MF 
observation in Ref.^^. The Majorana mode must be lo¬ 
calized both spatially and spectrally, i.e., the MF must 
be a well-defined zero energy excitation (in addition to 
being spatially localized at the wire ends) which becomes 
problematic at high temperatures where the zero mode 
hybridizes with the above-gap fermionic particle-hole ex¬ 
citations, thus completely losing its non-Abelian Majo¬ 
rana character. The energy localization property is inde¬ 
pendent of the spatial localization aspect, and both prop¬ 
erties must be satisfied for a system to have MF modes. 
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It has recently been showri^ that at such a high tem¬ 
perature the MF signature manifested in the tunneling 
experiment is diffuse over an energy range larger than the 
gap itself, and as such, the issue of MF localization be¬ 
comes moot because of the participation of other states. 
In the second part of our paper, we show that such broad 
and diffuse zero bias tunneling conductance peak could 
arise from subgap non-MF states which may generically 
be present in the system due to disorder. Thus, although 
the MF may indeed be strongly localized in the ferromag¬ 
netic nanowire system, tunneling experiments at temper¬ 
atures (and instrumental resolution) comparable to the 
gap energy cannot distinguish between MF features and 
ordinary (non-MF) subgap state features. Only future 
experiments carried out at much lower temperatures can 
therefore settle the question of what is being observed 
in the experiment of Ref.^. At much lower tempera¬ 
tures, however, the strongly localized nature of the MFs 
in the ferromagnetic nanowires will come into play in a 
dramatic fashion, leading to a strong zero bias conduc¬ 
tance peak in long wires and clearly split zero bias peaks 
in short wires, thus definitively establishing the existence 
(or not) of MFs in the ferromagnetic wire - superconduc¬ 
tor hybrid structure. On the other hand, if the physics of 
substrate-induced MF localization length suppression is 
not playing any role in the ferromagnetic adatom chains 
(as it is not in semiconductor nanowires) on supercon¬ 
ducting substrates, then at temperatures much lower 
than the induced gap, the subgap zero bias peak, if it 
is indeed arising from TS physics, should simply disap¬ 
pear completely since the MF localization length in such 
an unrenormalized situation would be much larger than 
the ferromagnetic adatom chain length one can fabricate 
on superconducting substrates at the present time. Thus, 
lowering (and improving) experimental temperature (in¬ 
strumental resolution) is the key to settling the question 
of whether MFs have indeed been observed in the exper¬ 
iments carried out in Ref.^. Since the state of the arts 
low-temperature STM experiments are routinely carried 
out at ~ 100 mK or below^ (a temperature regime acces¬ 
sible since the early 19903^, we urge future STM low 
temperature experiments (< 300 mK) in Fe chains on 
superconducting Pb substrates to settle the important 
question of the existence or not of MFs in this system. 


At this stage all one can say is that while Majorana 
modes do not show any basic inconsistencies with exper¬ 
iment, the broad features in the experimenlPS could also 
be consistent with accidental fermionic non-MF subgap 
states. In fact, the possibility that the observed sub¬ 
gap tunneling conductance structure in the experiment 
of Ref.l^ arises purely from a type of YSR bound states, 
rather than MF states, has been suggested recentlj@. 
In fact, the possibility that the observecP^^ broad and 
weak zero b ias feature arises from completely different 
physicd^^*^ nothing to do with Majorana fermions can¬ 
not be ruled out at this stage. 



FIG. 1. MF localization length (^) given by Eq. 0 as a 
function of dimensionless hopping strength (|t/Ap|) for vari¬ 
ous values of the chemical potential (|/r/Ap|). For |t| < |/r|/2 
the system is in the non-topological phase without Majorana 
fermions, and ^ is undefined in this regime. ^ diverges at 
|t| = |/r|/2, indicating the topological phase transition. For 
|t| ^ |/r|, C is well-approximated by |t/Ap|, which is the stan¬ 
dard coherence length formula for superconductors. 


II. MAJORANA DECAY LENGTH IN THE 
KITAEV CHAIN 


Let us first consider the prototypical and simplest 
model of a TS supporting Majorana end modes, the so- 
called Kitaev chain. This is a one-dimensional tight- 
binding model of spinless fermions with p-wave pairing, 
described by the Hamiltonian 

j 3 

A ^ ^ A h.c.) , (3) 

3 


where t is the hopping, fi the chemical potential, and Ap 
is the pairing potential. As shown by Kitae^P^, this model 
supports unpaired Majorana modes at its boundaries for 
|/x| < 2 |t|, with a MF localization length that is given by 




max 

A=±l 



Ay/-p 4A2 — p 
2{t+ |Ap|) 


-1 


( 4 ) 


We plot the localization length in Fig. as a function 
of the hopping amplitude for different values of p. Note 
that ^ is defined only for |t| > |p| /2, where the system is 
in the topologically non-trivial regime with a Majorana 
mode at each end. At \t\ = |p.| /2 the localization length 
diverges, indicating a topological phase transition into 
the topologically trivial regime at |t| < |p|/2. 

There are two special limits of interest. At the spe¬ 
cial point /r = 0 and |t| = |Ap|, the localization length 
vanishes and the MMorana is localized precisely at the 
end site of the chairP. We emphasize that in this fine- 
tuned case the localization of the Majorana is completely 
independent of the size of the energy gap, providing a 
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concrete example of a situation where a small gap could 
in principle also be associated with a small localization 
length. The other situation of interest is the physically 
realistic limit |f| ^ l^pl, where the bandwidth far ex¬ 
ceeds the superconducting gajP^. Here the localization 
length expressed to lowest order in Ap, reproduces the 
familiar form of a superconducting coherence length as 
discussed in the Introduction, 




t 


vf 

2i?gap ’ 


( 5 ) 


where Eg^p = |Ap|(l — and vp = 2|t|(l — 

are the spectral gap and Fermi velocity, re¬ 
spectively. Since |t| |Ap|, the localization length ^ ^ 1 

and the Majorana decay length is parametrically larger 
than the lattice constant (taken to be the unit of length 
here). On the other hand, it is clear that if for some 
reasons one can realize a Kitaev chain with |t| ~ |Ap|, 
as has been proposed for a quantum dot arrajEU, then 
the MF decay length is of order a few lattice sites only, 
qualitatively similar to the fine-tuned case. 


on symmetry considerations. Although spin itself is not 
a good quantum number in the presence of spin-orbit 
coupling, we can label the doubly-degenerate states by a 
pseudospin index ^ = ±1, which transforms like a spin 
under inversion and time-reversal. Assuming only a sin¬ 
gle band crosses the Fermi surface, the general expression 
for these states is 

\k,y= ^ [Bs,c;c,{k)\s,k,a) + 

( 6 ) 

where {k) and Bp^^„ (fe) are the coefficients of the 
s- and p-wave orbitals, respectively. Regarding the coef¬ 
ficients in Eq. § as 2 X 2 matrices in — cr space, one can 
derive a number of conditions. First, the normalization 
of the states in Eq. Q requires that 

so = Bl (fc) (fc) -f B; (fe) (fc). (7) 

Erom inversion and time-reversal symmetries we deduce 
that 

B, (fc) = Bs (-fc) = SyBl (-fc) Sy , (8a) 

Bp [k) = -Bp i-k) = SyB*p i-k) Sy . (8b) 


III. SUBSTRATE-INDUCED 
RENORMALIZATION OF THE TOPOLOGICAL 
WIRE 

We now turn to the physical realizatiorP^^^SMSHinilll 


In these equations Sp are Pauli matrices in the — a 
space. We additionally require that the pseudospin index 
behaves like a spin under mirror reflection in the planes 
perpendicular to the three Cartesian axes: 


of TS in a ferromagnetic nanowire in contact with a 

B, (fc)= 


= x,y,z, 

(9a) 

bulk s-wave superconductor. We note that the minimal 
model for the ferromagnetic wire proximity-coupled to a 
superconductoi*^l^l^AlMl^ as used in Ref.^^, is formally the 

Bp (fe) = 

f SpBp (Mpk) Sp 

p = x,y 

(9b) 


fl = z 


same as the corresponding semi conductor nanowire TS 
system introduced in Refs.®IISEIl.^^ith the only constraint 
being that the spin-splitting is arising from intrinsic ex¬ 
change splitting in the ferromagnetic system whereas it 
is induced as a Zeeman splitting by an externally ap¬ 
plied magnetic field in the semiconductor case. This has 
already been pointed out by Dumitrescu et aZ.I^. The 
presence of spin-orbit coupling and spin-splitting along 
with proximity-induced superconductivity enables us to 
avoid the fermion doubling theorem, leading to topologi¬ 
cal (effectively spinless p-wave) superconductivity in the 
wire. 


where Mp are reflection operators for the plane perpen¬ 
dicular to the /i-axis. 

In the pseudospin basis the bulk Green function of the 
superconductor is written as 


Gsc (w; k) = 


WTo -I- ^kTz + ^Scfx 


LO- 


_ ys _ a2 


( 10 ) 


sc 


where Tp are Pauli matrices in the particle-hole basis, 
is the dispersion in the superconductor, and Asc is the 
gap. The nanowire is placed on the (001) surface of the 
superconductor. The tunneling between the two systems 
is described by a superconductor-nanowire hopping term 


A. Self-Energy 

The normal state of the superconductor is character¬ 
ized by strong L • S spin-orbit coupling. Furthermore, two 
orbitals of different parity contribute to the states near 
the Fermi surface, which we label as s and p = p^, respec¬ 
tively. This is not essential for our theory (for example, 
we could also choose three distinct p-wave orbitals), but 
makes the following arguments more transparent; indeed, 
the details of the electronic structure of the supercon¬ 
ductor are not important for our analysis, which relies 


Btun — ^ ^ [^ncr {.^sSna T tpPna^ T h.C.j , (H) 

71(7 

where creates an electron with spin a at site n on 
the nanowire, and Sna- and Pna- are annihilation opera¬ 
tors of the orbitals of the superconductor. Only nearest- 
neighbor hopping is allowed. We remark here that a lo¬ 
cal breaking of inversion symmetry along z is required 
to generate non-zero tp which couples the nanowire sites 
to the p-orbital in the superconductor. Without loss of 
generality, we take ts and tp to be real. 
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The tunneling matrix T, which appears in the self¬ 
energy of the nanowire [Eq. 01. necessarily includes 
a transformation between the real-spin basis of the 
nanowire and the pseudospin basis of the superconduc¬ 
tor. It is written 


T = 


tsBT + 


i rT 

ip-Dp 


_ j. tdT _ j. pT 


( 12 ) 


With Gsc and T, the self-energy E can be readily com¬ 
puted as 


:(cc;r,rO= [ (k) Gsc iuJ-,k)T (k)^ 

\J ( Stt ) 


<pk 

= {1 — Z)ljJ + ZtTz + ZX ■ &Tz 


+ZAtx + ZA^*'^-&Tx, (13) 

where Z, t, A, A and are functions of w and r — r'. 
Their full expressions are given in the Appendix. 


B. Nanowire Hamiltonian 


We now consider the self-energy of the ferromagnetic 
nanowire in more detail. Assuming that the nanowire lies 
along the cc-axis, in the absence of the superconductor it 
can be modeled by the tight-binding Hamiltonian 


^int ^ h.C.) 


M ' ^ha-^ncT + B ^ ) 


^^ncr^niT ■ 


(14) 


where tint is the hopping intrinsic to the nanowire {not 
mediated by the superconductor), /r is the chemical po¬ 
tential, and B is the (spontaneous) ferromagnetic ex¬ 
change field (written out as an intrinsic magnetic field, 
rather than as an exchange splitting, in order to main¬ 
tain the explicit analogy to the semiconductor nanowire 
TS platforms where B is an extrinsic magnetic field). 

Including the self-energy due to the proximate super¬ 
conductor, the eigenenergies of the nanowire are given by 
the poles of the Green function 


Gwire (^) — 


1 


(15) 


^ — HBdG — S (w) UJ — H^S (uj) 
where TfudG is the BdG Hamiltonian of the bare 


nanowire Eq. (14), and in the second equality we have 


rearranged terms such that the effect of frequency renor¬ 
malization is captured by Z, and Hgg contains no terms 
proportional to lotq. Explicitly, 

77eff (^7 ^n) — ^ (<^m,n-t-l 4“ ^m,n—l) 

Z T ^ B6zn,n^z 


't~tm—n'^z 4“ ^m—r 

A*) 

^m—n ‘ X 


(TTz 


(16) 


The subscript m — n indicates that the quantities in 
Eq. (13) are evaluated at nanowire sites with relative 
coordinates r — r' = {xm — Xn) ^x- 

In general, the physics of the nanowire is ex¬ 
tracted from the Green function Gwire (w) including the 
frequency-dependent self-energy. Since we are interested 
only in the zero-energy Majorana mode and energy scales 
w <C Aq , however, we may take H^e (w = 0) as our effec¬ 
tive BdG Hamiltonian with no loss of generality. 


C. Effective Kitaev models 


To make further analytical progress we need to assume 
specific forms of Bg and Bp. We take 

Bg {k) = cos9so, (17a) 

Bp (^) — sin • (sy, Sxi ^.^o), (17b) 


where e^, = k/\k\. This choice is consistent with the 
symmetries of the pseudospin states Eq. and leads 
to an analytically tractable result which captures the es¬ 
sential physics we wish to explore. Other choices lead to 
qualitatively similar results. 

Using Eq. ( [T7| we calculate the full frequency- 
dependent forms of Z, t, A, A and which are given 
in the Appendix. At zero energy they take the relatively 
compact forms 


Z=l + 


tn>0 — 


Asc’ 
r cosna 




Z na 


'i>0 ' d* — 1(7yQ 


sin 20 


r cos na 4- na sin na 


A„ = 


r sinna 

7 

Z na 


-n/i 


Ai>o ■ ^ = *o-yffsin20 


r sin na — na cos na 
Z n^a? 


(18a) 

(18b) 

-"(f8c) 

(18d) 

^/«(,18e) 


where we have introduced dimensionless variables 


T = Kv {t1 cos^ ^ 4- sin^ 0) , (19) 

g = TTvtgtp, a = /cFfliat, and ^ = C/aiat, in which v 
and ^ are respectively the Fermi-level density of states 
and the coherence length of the superconductor, and aiat 
is the lattice constant of the tight-binding model of the 
nanowire [Eq. fill- 

in the limit of large exchange field (H 3> |A„|, | A„|), 
the effects of the spin-orbit coupling and s-wave pair¬ 
ing terms are suppressed. As shown in Appendix [B| the 
Hamiltonian ( [I^ thus reduces to two copies of the Kitaev 
model, albeit with long-range hopping. To make connec¬ 
tions with Sec. HD we first ignore the long-ranged part of 
the self-energy which is beyond nearest neighbors, yield¬ 
ing (second quantized) effective Hamiltonians for spin-up 












(+) and spin-down (—) species 

Hq ^ = — (Z ^^int + ii) + h.c.^ 

j 

-Z-^ a]aj 

3 

-|-AJ ^ 'y ^ (a^cij+i -|- h.c.). (20) 

3 


The induced hopping integral and pairing potential are 
given by 




r COSO n/f 

Z a 



. r sin a — a cos a 

0‘"'-‘‘z -- 




( 21 ) 

( 22 ) 


From Eq. the localization length for the Majorana 
zero modes in these Hamiltonians (valid in their topolog¬ 
ical phase) is thus 


^wire 


Zi ^tint + tl 


r,(±) 


2E, 


(±) 


(23) 


where the renormalized Fermi velocity and excitation gap 
are 


r,(±) 


— 2 (Z ^tint + il) A/l ~ 




4 (tint + ^tl) 


(24) 




i^^T BY 


4 (tint + Zti) 


(25) 


Note that the localization length is the same for the spin- 
up and -down sectors. From Eq. (23) we observe that 
if one ignores the renormalization of the Fermi veloc¬ 
ity and uses instead its intrinsic value, = 2tint[l — 
(/i =F bY /4t?jt]^/^, to estimate the localization length as 


5 = u^V2^'gap\ the result would overestimate the true 

value by a factor of > 1- If fh® coupling 

between the wire and the superconductor is weak (i.e. 
F <C Age, tint), the velocity is only weakly renormalized 
and However, when F is comparable to Agq 

or even tint, the discrepancy between the renormalized 
and the bare Fermi velocity is huge. For large enough F 
and hence Z, the coherence length could be close to zero 
even though the induced triplet gap is small. Whether 
or not this strong velocity renormalization, leading to 
sharply-localized MFs in the TS nanowire, is present in 
the experiment of Ref.^ can only be determined empir¬ 
ically since the microscopic details about F are simply 
not known in the experimental system. What is clear, 
however, is that there is a well-defined physical mecha¬ 
nism, namely, a very strong tunnel-coupling between the 
superconductor and the nanowire, which would lead to 
a strong renormalization of the effective Fermi velocity 
and a concomitant suppression of the MF localization 
length in the nanowire even if the induced topological 


gap is small. We note that the existence of the strong 
renormalization effect has already been invoked for the 
Fe/Pb system by Peng et al. using a helical magnetic 
chain model for the nanowir^^. 


D. Effects of non-local hopping and pairing 


As we mentioned in the introduction, the substrate- 
induced enhancement of the Majorana localization length 
is accompanied by a power-law decay of the MF^^. This 
power-law decay of the MFs, if large, limits the validity 
of the enhanced exponential localization. To understand 
and estimate this effect we write the Hamiltonian in the 
large tunneling limit as 


H = Ho+ SH 


(26) 


where Hq is given in Eq. (201 and SH contains the hop¬ 


ping and pairing terms in Eq. (16) involving sites sepa¬ 


rated by two or more lattice spacings. Let ipo denote the 
zero-energy Majorana mode that is localized at the end 
of the wire with a localization length given by Eq. (231. 
With the non-local perturbation the state acquires 

a correction -ipo ^ Yo = Yo + SYo, where 


6Yo = -^V5H\Yo) , 
Bo 


(27) 


where V = 1—IV'o) (f/’ol- We can now qualitatively see the 
localization behavior of Yo including the long-range self¬ 
energy correction. We begin by re-writing the equation 
for 6Yo in real space with coordinate x as 


8Yo{x) 


dx'dxiGoix, x')8H{x', xi)Yo{xi ), (28) 


where Go{x,x’) « HY^V is the projected Green func¬ 
tion corresponding to the Hamiltonian Hq at zero energy. 
Eor a gapped system the Green function Go{x,x') oc 
g-k-a; |/?wiro unperturbed wavefunction Yo 

g-a:/Cwire figth localized on a length scale ,^wire- On 
the other hand, the perturbation 6H has a power-law tail 
so that 5H{x', xi) oc J[kpW — xi\) where ^ is 

the coherence length of the bulk superconductor. Assum¬ 
ing that ^wire “C C) we see that the integral in Eq. is 
dominated by |xi|, |a; — x'| < Cwire so that for x ^ ^wire 
the perturbed wave-function is written as 

= V’o + ^ + P- -, (29) 

kpx 


where a is a normalization constant and /? is a param¬ 
eter determined from the perturbation theory. Strictly 
speaking, the localization length of Yo is max{^^ire,C}- 
However, the factor of kp in the denominator of the sec¬ 
ond term in Eq. (291 means that this component of the 
wavefunction is effectively localized on a scale of ^ Y 
This recalls YSR bound states around magnetic impuri¬ 
ties in conventional superconductoil^, where the bound 
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state wavefunction has a similar spatial profile in the 
superconductor. Indeed, experiments show that such 
bound states are localized on a scale much smaller than 


The implication of Eq. (29) is that since the dominant 


localization length of the non-local part is small, the non¬ 
local term 5H can be safely ignored: the experimentally 
measured localization length will still be ^wire even with 
longer-range hopping in Eq. |16[ We note, however, that 
independent of whether a ^ (the perturbative regime) 
or P ^ a (the YSR regime) in Eq. (29), the resultant Ma- 


jorana wavefunction is strongly localized at the wire end 
{x = 0) with a localization length of ^wire or , both 
of which are much smaller than the bare ME localization 
length ^ expected without the substrate renormalization 
effect (provided, of course, one is in the strong tunnel 
coupling regime). In Ref.“ Dumitrescu et al. recently 
took into account the second term in Eq. (291 as causing 


the suppressed ME localization in ferromagnetic chain 
TS systems whereas Peng et al!^ mostly considered the 
first term in discussing ME localization in helical mag¬ 
netic chains. In principle, both terms could be important, 
but their qualitative effects are similar, both leading to 
a strongly suppressed ME localization in the nanowire 
in the strong tunnel-coupled regime. Importantly, the 
power law decay of the second term in Eq. (29) may have 


negative implications for the non-Abelian braiding exper¬ 
iment^^, severely limiting the usefulness of the resultant 
MFs in carrying out topological quantum computation 
although for practical purposes the MFs appear strongly 
spatially localized at the wire ends. 


E. Relating quasi-lD models to ID models 


We have established above that as long as the nanowire 
is strongly tunnel-coupled to the superconductor (so that 
the condition E > tint,Asc applies), the ME localiza¬ 
tion length would be strongly suppressed compared with 
the standard bare coherence length formula due to the 
Fermi velocity renormalization caused by the substrate. 
This renormalization effect appears to be independent 
of the nature of the nanowire and, therefore, should 
affect both ferromagnetic nanowires and semiconductor 
nanowires equally (as long as the tunnel coupling defined 
by Eq. (191 is large). We now show that this is not the 
case, and there is good reason to believe that the ferro¬ 
magnetic chain system of Ref.^ is much more strongly 
renormalized by the substrate than the semiconductor 
nanowire system^^SHlol 


While it is common to assume a strictly ID limit in 
modeling both the semiconductor nanowire and ferro¬ 
magnetic chain systems, a more realistic model would 
treat them as quasi-ID and as a result the parame¬ 
ters ts and tp in the ID model [for example, defining 
r in Eq. (EU)] are really effective parameters that have 
a strong dependence on the radius of the nanowire in 
a quasi-lD geometry. Since we are interested in un¬ 


derstanding the scaling behavior with nanowire radius, 
we will assume a simple model of a three-dimensional 
(3D) cylindrical lattice nanowire (with the wire trans¬ 
verse cross-sectional width being much smaller than the 
wire length). The 3D (i.e. quasi-lD) wavefunctions and 
the strictly ID wavefunctions '(/;™®’^'^(a::) are related by 
a transverse wavefunction factor as 


/ wire,3D / / \ 


'yf7r(R -t- u)(/cm ) 

Xgim0,^wire.lD 


Jjr 


P 


(30) 


where km,j is the zero of the Bessel function so 


that the wavefunction satisfies '0. 


wire,3D 


{R + a,x,(j)) = 0 


and a is the lattice constant of the wire. In deriv¬ 
ing this expression we have used the continuum limit 
kmja/R ^ 1 so that the lattice wavefunction can be 
approximated as a times the continuum wavefunction 
of a cylinder. Requiring that the lattice wavefunction 
vanishes at the edge of the cylinder corresponds to the 
boundary condition that the continuum wavefunction 
vanishes a distance a outside the wire. 

The ID hopping matrix elements enter the for¬ 
malism in Sec. |III C] through the parameter T defined in 
Eq. 19 To simplify our analysis we split T = Tg cos 9'^ -I- 


Tp sin0^ where Tg = TTvtl and Tp = 'Kvt^. The matrix el¬ 
ements tg^p are proportional to the absolute value square 
of the 3D tunnelings ta^p and the quasi-lD wavefunction 
{p,x,(j)) at the surface of the wire. In addition, 
for the purpose of our estimate, we will make a simplify¬ 
ing assumption that the Green function of the substrate 
(which determines the density of states v = Im{G}) is 
local, i.e. kpa ^ 1. With these assumptions (which may 
be relaxed without qualitatively changing the results), 
the generalized form for Tg^p is written 


r^.p = TrnFa,p 2:, 0) p 

surf 

« J {R, X, 0) p, (31) 

where the sum in the first line is over all lattice sites at 
the surface of the wire. For R ^ a the wavefunction at 
the surface is approximately given by 

, (32) 

and so we find 


^ S,p '^‘‘s,p ^3 ■ 


(33) 


The hopping fg p can be parametrized by a dimensionless 
parameter (^s,p and written as ts,p = Cs,pEf. Since the 
hopping between the wire and the SC is expected to be 
of the same order of magnitude as the bare hopping in 
the substrate, we expect (^s,p to be a parameter of order 
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1. while the density of states is ty = kpa^ Using 
this parametrization we estimate as 

(34) 

where k = (g pk^ j{kFaY/AF‘^. 

We are now in a position to compare F for the semi¬ 
conductor nanowire and the ferromagnetic chain. Quali¬ 
tatively speaking, R is much smaller for a ferromagnetic 
chain of atoms (e. g^_^Ref.l^ than for the semiconductor 
nanowire (e.g. RefP^, and so F is accordingly expected 
to be much larger in the former compared to the latter 
due to the dependence F oc R~^. Quantitatively, as¬ 
suming i? ^ a is of order 0.5 nm for the ferromagnetic 
Fe chain in Ref.l^, we estimate F ^ Ep leV (assum¬ 
ing K ~ 1). On the other hand, for the semiconduc¬ 
tor nanowire with much larger mode confinement radius 
R ~ 20nm, the same values for the other parameters 
yields a self-energy parameter F on the order O.OlmeV. 
This huge difference in F between the semiconductor 
nanowires used iiPl 

and the ferromagnetic chains used 
iipa gives a natural explanation for why the MF might 
be strongly localized (delocalized) in the ferromagnetic 
(semiconductor) nanowires even if both systems mani¬ 
fest the same induced superconducting gap (~ 100/xeV). 
This difference arises, keeping all the other parameters 
similar, from the difference in the transverse quantiza¬ 
tion size in the two systems. The wire radius ratio of 
roughly a factor of 40 between them can in principle lead 
to a localization length difference by as large as a factor 
of 40^ = 64000! In reality, this is an overestimate, since 
the bare Fermi velocity in the semiconductor is typically 
a factor of 100 or so smaller than that in the ferromag¬ 
netic metallic chain and also the tunneling factor (^g p 
in the semiconductor can easily be smaller by an order 
of magnitude. The combination of these factors may re¬ 
duce the factor from 64000 to 64000/100 ~ 640 difference 
in the MF localization length between the semiconduc¬ 
tor nanowirdi^lt^o] ferromagnetic wir^^ systems 

even if both have exactly the same induced supercon¬ 
ducting gap (~ 0.1 meV). This factor of ~ 500 difference 
is in quantitative agreement with the conclusion of Ref. 
21 where the MF localization length is inferred to be 
< 1 nm whereas in the semiconductor nanowire case the 
MF localization length is the same as the bare coherence 
length in the nanowire (^ 100nm). Thus, the differ¬ 
ence between MF localization in the two systems arises 
entirely from the difference in the nanowire transverse 
confinement radius in semiconductors versus metals. 

In the next section we show that finite temperature 
effects in the experiment oP^^ would make the Majorana 
zero-mode signature weak and diffuse in the tunneling 
conductance measurement (in spite of strong MF local¬ 
ization) because of strong thermal broadening of both 
the MF mode and above-gap fermionic excitations in the 
system. Thus, the high-temperature MF signature (at 
temperatures comparable to the gap energy) is qualita¬ 
tively similar in the tunneling spectroscopy as that of 


any generic non-MF subgap excitation even when these 
fermionic subgap excitations are not at zero energy. Only 
experiments at temperatures much lower than the in¬ 
duced gap energy can distinguish MF versus ordinary 
fermionic subgap states in the tunneling spectra indepen¬ 
dent of the localization properties of the Majorana zero 
modes. At temperatures much lower than the gap energy, 
however, the strongly suppressed MF localization length 
in the ferromagnetic nanowire becomes an extremely im¬ 
portant physical quantity since the very short MF local¬ 
ization length may now allow well-defined (rather than 
strongly hybridized) zero-energy MF modes to exist in 
rather short magnetic chains used in Ref.^, which would 
not be possible if the MF localization length is given by 
the bare formula. Whether this physics is operational or 
not can only be determined empirically by carrying out 
STM measurements at temperatures much lower than the 
gap energy. 


IV. IMPURITY-INDUCED SUBGAP STATES 

Split Majoranas that may result from a wire of length 
comparable to the localization length are essentially iden¬ 
tical to conventional Andreev states that are of non- 
topological origin. In this section we examine tight- 
binding models for the tunneling conductance from an 
STM tip into either a p-wave or an s-wave nanowire 
with subgap states of nontopological origin. In the p- 
wave case, we induce subgap Andreev states by includ¬ 
ing non-magnetic impurities. In the s-wave case, on the 
other hand, YSBP^ subgap states are created by mag¬ 
netic impurities. The calculated finite temperature tun¬ 
neling conductance results (Fig. ^ due to these fermionic 
subgap states are compared with the corresponding MF- 
induced tunneling spectra in a TS (Fig. |^. We demon¬ 
strate that, in general, the two sets of results are almost 
impossible to distinguish unless the experimental tem¬ 
perature and energy resolution are much smaller than the 
superconducting gap energy. We mention here that the 
temperature dependence of the MF tunneling conduc¬ 
tance spectra for the ferromagnetic nanowire system has 
already been calculated in great details by Dumitrescu 
et al^, who have established that high-temperature 
ZBCP is spectrally spread over the whole energy gap as 
a very weak and very broad feature making the inter¬ 
pretation of the data in Ref.^ problematic. Our results 
for the temperature dependence of MF induced ZBCP in 
the ferromagnetic nanowires agree completely with the 
results presented in Ref.^, but what is new in our cur¬ 
rent work is showing that non-MF subgap states may also 
lead to tunneling conductance features which are indis¬ 
tinguishable from the corresponding MF features in high 
temperature experiments. This comparison between MF 
versus non-MF conductance features is the new ingredi¬ 
ent in our results. 
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A. Model for a p-wave nanowire 

The Hamiltonian describing tunneling into the p-wave 
nanowire is 


H — Hnj + Ulead -^tun 

(35a) 

Hw = (a]aj+i -bh.c.^ 

j 

j 

+ Ap ^ ^ + h.C.) 

j 

(35b) 

-^tip — ^ ^ 
k 

(35c) 

Htnn = tpclai + h.C. 

(35d) 


where H^, i?tip and i?tun are respectively the Hamilto¬ 
nians for the nanowire, the STM tip, and the tunneling 
from the tip to the nearest site on the nanowire (desig¬ 
nated as site 1). The annihilation operator for site j in 
the nanowire is aj, while the annihilation operator for 
state k in the STM tip is Cfc. The p-wave spectral gap of 
the nanowire is A = Ap(l — fj?. Non-magnetic 
impurities are introduced by giving the chemical poten¬ 
tial a site-dependence. We ignore spin as we are only 
interested in non-magnetic impurities. 

The zero-temperature conductance is computed from 
the Green function Q (x,x',uj) of the nanowire bjUSl 


Gt=o (w) = ^ [27r7pi (w) - 7^ \gl^ (w)| -b 7 ^ |/i\ (w)|] 

(36) 

where g and / are respectively the normal and anoma¬ 
lous parts of Q, and the subscripts “ 11 ” indicate that 
both position arguments of the Green function are at 
the nanowire site in contact with the STM tip. 7 = 
2TTvt\ is the broadening due to the STM tip and pi = 
— (l/ 7 r)Im 5 ii is the local density of states at the point of 
contact with the tip. Numerically, the Green function is 
given by 


Q {x, x', oj) 


uj + iS — Hyj (x, x') -b i-F (x, x') 


(37) 

where T (x, x') = ^5x,i5x',i is the broadening induced by 
the contact to the STM tip. The broadening term iS 
mimics the energy resolution of the setup, with a lower 
5 implying a higher resolution. We note that the experi¬ 
ment of Ref.^ has a very poor energy resolution, compa¬ 
rable to the induced energy gap in the nanowire, making 
the broadening parameter 5 an important aspect of the 
experimental analysis. 

The finite-temperature conductance can be obtained 
from Gt=o using 


Gt (w) = J deGr^o (s) d^np (w — e) 

= ^/ deGT=o (e)sech^ 


(38) 


where np (oj) is the Fermi distribution function. The 
thermal broadening of the zero-temperature conductance 
is also important to understanding the high-temperature 
STM experiments carried out in Ref.l^ as already em¬ 
phasized in Ref.lSI. 


B. Model for an s-wave nanowire 

The Hamiltonian for tunneling into the s-wave 
nanowire is 


H — Hyj -b Fflead + Htun (39a) 

ddw = ~t 'y ^ ^Oj^ajpi^cr + h.c.^ — /i ^ 

jo- j<7 

+ (39b) 

j 

Fdtip — ^ (39c) 

ka 


ddtun — 'y ( “b h.C. 

(7 


(39d) 


We again adopt the convention that the nanowire site 
in contact with the STM tip is denoted as site j = 1. 
The definition of the annihilation operators is general¬ 
ized to include the spin degrees of freedom, which must 
be accounted for in this case. In contrast to the p- 
wave nanowire, here we take a uniform chemical poten¬ 
tial p, but allow for the possibility of magnetic impurities 
through the site-dependent Zeeman field Bj. The spec¬ 
tral gap is A = Ag. 

The analysis of the tunneling proceeds similarly to that 
for the p-wave nanowire above. The expression for the 
zero-temperature conductance is, however, slightly mod¬ 
ified to include the spin degrees of freedom 


Gt=o (w) = 


27r7^Pi,T (w) 

a 


(40) 


where g and / are the normal and anomalous Green func¬ 
tions, respectively, and the subscripts indicates the site 
and spin indices. 


C. Effects of high temperature and low resolution 
on the tunneling conductance 

Fig-i shows the differential conductance obtained by 
tunneling into nontopological subgap states in an infinite 
nanowire, induced by either magnetic impurities (in the 
s-wave nanowire) or non-magnetic impurities (in the p- 
wave nanowire). At high resolution and low temperature 
(Fig. [2}r ,d,g,j) the exact non-zero energy of the subgap 
state can be extracted from the tunneling spectra. With 
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w/A w/A w/A 



FIG. 2. (a-c) The conductance measured by placing the STM tip at a site with a magnetic impurity in an infinite s-wave 
nanowire with t = — lOA and 7 = O.OOIA. The magnetic impurity induces a YSR state with subgap energy 0.5A. We show 

the conductance at high resolution and low temperature [panel (a)], high resolution and high temperature [panel (b)[, and low 
resolution [panel (c)]. (d-f) Similar to Panels (a-c) but with smaller subgap energy 0.2A. (g-i) Similar to Panels (a-c) but at 
the site of a nonmagnetic impurity in an infinite p-wave nanowire with f = 2p = lOA. The impurity induces an Andreev state 
with energy 0.5A. (j-1) Similar to Panels (g-i) but with smaller subgap energy 0.2A. 


either high temperatures (Fig. ,h,k) or low resolu¬ 
tions (Fig.[^,f,i,l), however, the peaks broaden and could 
appear to arise from a broadened zero-energy peak. In 
fact, Fig. shows that with high temperature or low res¬ 
olution (or both), the tunneling conductance results in 
the nanowires typically manifest broad features consis¬ 
tent with a zero-energy peak as long as the subgap states 
are located at \E\ < 0.5A. Thus, the observation of broad 


zero-bias conductance features should not be associated 
as evidence for the existence of precisely zero-energy MFs 
since this is also consistent with fermionic subgap states 
such as ordinary YSR or Andreev bound states. Indeed, 
the possibility that the experiment of Ref.^ is actually 
observing a YSR state feature instead of a MF state has 
recently been suggested in the literatur^^. 

We show in Fig. [^a)-(d) the conductance of a genuine 
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FIG. 3. The differential conductance for an STM tunneling 
measurement on a p-wave nanowire with t = 2/i = lOA and 
7 = O.OOIA. In panels (a)-(d), the STM tip is placed at the 
end of the wire, where a zero-energy Majorana mode is found. 
The results are computed under the conditions of: (a) high 
resolution and low temperature, (b) high resolution and high 
temperature, (c) low resolution and high temperature, and 
(d) low resolution and high temperature. In panels (e)-(f), the 
STM tip is placed at a site with an impurity, which induces 
a subgap state with energy 0.3A, and the results shown cor¬ 
respond to (e) high resolution and (f) low resolution, respec¬ 
tively. The high-temperature results in (b,d,f) are very sim¬ 
ilar, indicating the difficulty in distinguishing between states 
at zero energy and states at non-zero but small energy in such 
circumstances; on the other hand, low-temperature conduc¬ 
tance curves with Majoranas in (a,c) can be clearly distin¬ 
guished from the curves in (e). 


Majorana mode at the end of a clean p-wave nanowire, 
with no other subgap states being present in the system. 
For comparison, in Fig.[^e)-(f) we show the conductance 
of a subgap non-MF state induced by a non-magnetic 
impurity in an otherwise uniform p-wave nanowire. We 
observe that for low resolutions or high temperatures, 
the conductance spectra for zero-energy Majorana modes 
cannot be effectively distinguished from those for non- 
MF subgap modes. In particular, the high temperature 
(and/or poor energy resolution) plots in Figs. I^e) and 
(f) are indistinguishable from the corresponding plots 
in Figs, [^b)-(d), thus confirming that high-temperature 
STM experiments, as carried out in Ref.“ cannot really 
confirm the existence of Majorana modes. Thus, only 
future experiments at lower temperatures and employing 
better instrumental resolution would be able to conclu¬ 


sively determine the existence or not of MFs in the hybrid 
Fe nanowire-Pb superconductor system recently studied 
in Ref.1211. 

Our work (as well as that presented in Ref.^^ shows 
that any subgap states, whether MF or not (and whether 
precisely at zero-energy or not), would manifest very sim¬ 
ilar STM conductance spectra at high experimental tem¬ 
peratures, and thus the observation of broad, weak, and 
spectrally diffuse zero bias conductance “peak” cannot 
be taken as being synonymous with the existence of MFs 
in the system, as has been done rather uncritically in 

Ref.l^. 


V. CONCLUSION 

In summary, we have established that the nanowire 
on superconductor hybrid systems can potentially have 
very short Majorana localization length even when the 
induced topological superconducting gap is very small in 
the nanowire by virtue of the substrate induced strong 
renormalization of the effective nanowire parameters (e.g. 
the Fermi velocity, the gap, etc.) because of strong 
frequency dependence of the relevant self-energy func¬ 
tion determining the proximity-induced pair potential 
in the nanowire. We have shown that this renormal¬ 
ization goes as R~^ where R is the effective nanowire 
confinement size in the transverse direction determin¬ 
ing how one-dimensional the system really is (with R 
going to zero limit being the true ID nanowire limit). 
This provides an explanation for why the Majorana lo¬ 
calization length could be very small (large) in metal¬ 
lic (semiconducting) nanowires on superconductors since 
R ^ 0.5(20) nm in the two systems leads to a large dif¬ 
ference in the renormalization effect induced by the sub¬ 
strate. The substrate-induced suppression of the Majo¬ 
rana localization length may have implications for recent 
effort^^ to observe localized Majorana modes in fairly 
short (;< 50nm) ferromagnetic Fe chains on supercon¬ 
ducting Pb substrates usiM STM spectroscopy, provid¬ 
ing a possible explanatiorP^ for how the Majorana mode 
may be spatially highly localized on a sub-nm length scale 
near the ends of the Fe adatom chain in spite of a very 
small induced superconducting gap. 

But, the definitive observation of spatially localized 
Majorana fermions (rather than just spatially localized 
ordinary fermionic subgap states) requires precise energy 
localization (exactly at mid-gap or zero energy) in addi¬ 
tion to strong spatial localization. Such an “energy local¬ 
ization” necessitates an experimental temperature much 
smaller than the topological superconducting gap. This 
condition of “energy localization” is absent in the experi¬ 
ment of Ref.since the temperature is larger than the es¬ 
timated induced superconducting gap. Thus, the spectral 
weight of the observed “zero bias” conductance feature is 
spectrally spread over the whole induced gap, making it 
difficult to conclude about possible existence of Majorana 
fermions in spite of the observed spatial localization. In 
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addition, we have shown through explicit numerical sim¬ 
ulations that high-temperature and low-resolution tun¬ 
neling conductance measurements cannot distinguish be¬ 
tween Majorana modes and ordinary fermionic subgap 
states as both manifest broad and weak zero-bias con¬ 
ductance features. The distinction between these two 
situations necessitates experiments at temperatures (and 
resolutions) well below the induced superconducting gap 
energy in the nanowire. We mention in this context that 
the state of the arts STM experiments on superconduc¬ 
tors are routi nely carried out at temperatures of 100 mK 
or belov02ll31^ and future experiments with much lower 
temperatures (and much better resolutions) than used 
in Ref.^ can settle the question of whether Majorana 
fermions have indeed been observed or not in ferromag¬ 
netic chains on superconducting substrates. The mini¬ 
mal existence proof of Majorana fermions necessitates the 
demonstration of both spatial (at wire-ends) and spectral 
(at zero-energy) localization of the observed excitation. 

Lowering temperature (and/or enhancing the topolog¬ 
ical gap) as well as improving instrumental resolution 
should lead to the MF-induced ZBCP becoming sharper 
and stronger [as in Fig. [^a)-(d)] in longer ferromagnetic 
chains whereas in shorter chains, the ZBCP should split 
because of the hybridization between the two Majorana 
end mode^^mm (The currently observed zero bias peak 
in the ferromagnetic chain^^ is only a very weak frac¬ 
tion (10“^) of the quantized value of 2e^//i predicted for 
the Majorana zero modes, and lowering of temperature 
should enhance its strengtlP^.) On the other hand, if 
the observed sub-gap STM conductance features are aris¬ 
ing from impurity-induced non-MF subgap bound states 
[as in Fig. or Fig. [^e)-(f)], then lowering tempera¬ 
ture (and/or increasing the topological gap) and enhanc¬ 
ing resolution should clearly show that these accidental 
subgap states are non-topological and therefore the re¬ 


sultant subgap conductance features are not spectrally 
located exactly at zero energy. There should not be any 
energy splitting of such non-MF subgap peaks in shorter 
chains, also distinguishing them from possible MF peaks. 
Current experiments^ do not actually find any clear ev¬ 
idence for a topological gap in the ferromagnetic chain 
with the background subgap conductance being typically 
more than 50% of the normal state conductance, and any 
definitive claim for the observation of Majorana fermions 
must necessarily be coupled with the observation of a 
reasonably well-defined topological gap with background 
subgap conductance being suppressed by at least one or 
two orders of magnitude from its current value. In this 
context, it may be useful to remember that any fermionic 
subgap bound state in the adatom (Fe)-superconductor 
(Pb) system will most likely also be localized near the 
ends of the Fe chain since this is where the impurity po¬ 
tential is the strongest, as was explicitly shown recentlj®^ 
for possible YSR states in this system. Indeed a direct 
numerical diagonalization using model band structures 
of the Fe/Pb experimental system produces a very large 
complex of states in the superconducting gajP^, most of 
which have nothing topological about them (but all of 
which are likely to contribute to the subgap tunneling 
conductance), and the complexity of the system there¬ 
fore makes any straightforward interpretation of the STM 
data very difficult. 
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Appendix A: Derivation of Effective Hamiltonian 


The various terms in the expression of the self-energy [Eq. @] are 
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where Mi = +tpBpB* and M 2 
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If we take the specified form of Bg and Bp [Eq. (17)], 
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Further, restricting r only at the lattice positions of the nanowire lying on the x-axis, we can evaluate the terms 
exactly as 
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The full effective BdG Hamiltonian of the nanowire is then given by substituting these expressions to Eq. (16). 
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Appendix B: Hamiltonian in the strong spin 
polarization limit 

To obtain the Hamiltonian in the limit of strong spin 
polarization, we start by writing Eq. (16) in momentum 
space 

HeS = ^ |^(efc - palak + Xto-layak + Saj^CT^Ofe 

k 

+Afc^ (afcT“-fcT + + ^-c-) (Bl) 

+Afc (afefO-fe^ + h.c.)], 

where a\ = efe is the Fourier transform of 

the hopping terms in Eq. (16), B = Z~^B is the renor¬ 


malized exchange field, while Afe, Afc, and A^*^ are the 


proximity-induced spin-orbit coupling, singlet pairing po¬ 
tential, and triplet pairing potential, respectively. Di- 
agonlizing the first three terms in Eq. (Bl) gives us the 


“normal state” dispersion ^±{k) = Ck — fJ, i y B"^ + \\ 
with eigenvectors 
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■ (B2) 


Eollowing the standard approaclF^, we rewrite the 
Hamiltonian in this basis with the transformation 


Ofc = (j)+ (fe) C+ (fe) -I- 4>- (fe) c_ (fe) , 


(B3) 


where c± annihilates states in the corresponding normal 
state bands. Eq. (Bl) then becomes 
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The pairing terms on the first line correspond to the 
proximity-induced triplet gap. On the second line, we 
have additional triplet pairing terms induced from the 
interplay of the proximity-induced spin-orbit coupling Afc 
and singlet gap A^. In the strongly spin-polarized limit 
where B ^ |Afc|, |Afe|, however, these terms are negli¬ 
gible. The singlet pairing term, also on the second line. 


is not suppressed at the Hamiltonian level. However, in 
the strongly spin polarized limit the band basis approx¬ 
imately coincides with the spin basis, and so the singlet 
pairing potential does not open a gap due to the huge mo¬ 
menta mismatch between the two bands, and may thus 
be ignored in the effective Hamiltonian. 
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